Unraveling the genetic architecture of congenital vertebral malformation with reference to the developing spine.

Congenital vertebral malformation, affecting 0.13-0.50 per 1000 live births, has an immense locus heterogeneity and complex genetic architecture. In this study, we analyze exome/genome sequencing data from 873 probands with congenital vertebral malformation and 3794 control individuals. Clinical interpretation identifies Mendelian etiologies in 12.0% of the probands and reveals a muscle-related disease mechanism. Gene-based burden test of ultra-rare variants identifies risk genes with large effect sizes (ITPR2, TBX6, TPO, H6PD, and SEC24B). To further investigate the biological relevance of the genetic association signals, we perform single-nucleus RNAseq on human embryonic spines. The burden test signals are enriched in the notochord at early developmental stages and myoblast/myocytes at late stages, highlighting their critical roles in the developing spine. Our work provides insights into the developmental biology of the human spine and the pathogenesis of spine malformation.

Data Mining Suggests That CXCL14 Gene Silencing in Colon Cancer Is Due to Promoter Methylation.

CXCL14 is one of the most evolutionarily conserved members of the chemokine family and is constitutionally expressed in multiple organs, suggesting that it is involved in the homeostasis maintenance of the system. CXCL14 is highly expressed in colon epithelial cells and shows obvious gene silencing in clinical colon cancer samples, suggesting that its silencing is related to the immune escape of cancer cells. In this paper, we analyzed the expression profiles of multiple human clinical colon cancer datasets and mouse colon cancer models to reveal the variation trend of CXCL14 expression during colitis, colon polyps, primary colon cancer, and liver metastases. The relationship between CXCL14 gene silencing and promoter hypermethylation was revealed through the colorectal carcinoma methylation database. The results suggest that CXCL14 is a tumor suppressor gene in colorectal carcinoma which is activated first and then silenced during the process of tumor occurrence and deterioration. Promoter hypermethylation is the main cause of CXCL14 silencing. The methylation level of CXCL14 is correlated with the anatomic site of tumor occurrence, positively correlated with patient age, and associated with prognosis. Reversing the hypermethylation of CXCL14 may be an epigenetic therapy for colon cancer.

Variants in the SOX9 transactivation middle domain induce axial skeleton dysplasia and scoliosis.

SOX9 is an essential transcriptional regulator of cartilage development and homeostasis. In humans, dysregulation of SOX9 is associated with a wide spectrum of skeletal disorders, including campomelic and acampomelic dysplasia, and scoliosis. The mechanism of how SOX9 variants contribute to the spectrum of axial skeletal disorders is not well understood. Here, we report four novel pathogenic variants of SOX9 identified in a large cohort of patients with congenital vertebral malformations. Three of these heterozygous variants are in the HMG and DIM domains, and for the first time, we report a pathogenic variant within the transactivation middle (TAM) domain of SOX9 . Probands with these variants exhibit variable skeletal dysplasia, ranging from isolated vertebral malformation to acampomelic dysplasia. We also generated a Sox9 hypomorphic mutant mouse model bearing a microdeletion within the TAM domain ( Sox9 Asp272del ). We demonstrated that disturbance of the TAM domain with missense mutation or microdeletion results in reduced protein stability but does not affect the transcriptional activity of SOX9. Homozygous Sox9 Asp272del mice exhibited axial skeletal dysplasia including kinked tails, ribcage anomalies, and scoliosis, recapitulating phenotypes observed in human, while heterozygous mutants display a milder phenotype. Analysis of primary chondrocytes and the intervertebral discs in Sox9 Asp272del mutant mice revealed dysregulation of a panel of genes with major contributions of the extracellular matrix, angiogenesis, and ossification-related processes. In summary, our work identified the first pathologic variant of SOX9 within the TAM domain and demonstrated that this variant is associated with reduced SOX9 protein stability. Our finding suggests that reduced SOX9 stability caused by variants in the TAM domain may be responsible for the milder forms of axial skeleton dysplasia in humans.

Dissection of mendelian predisposition and complex genetic architecture of craniovertebral junction malformation.

The craniovertebral junction (CVJ) is an anatomically complex region of the axial skeleton that provides protection of the brainstem and the upper cervical spinal cord. Structural malformation of the CVJ gives rise to life-threatening neurological deficits, such as quadriplegia and dyspnea. Unfortunately, genetic studies on human subjects with CVJ malformation are limited and the pathogenesis remains largely elusive. In this study, we recruited 93 individuals with CVJ malformation and performed exome sequencing. Manual interpretation of the data identified three pathogenic variants in genes associated with Mendelian diseases, including CSNK2A1, MSX2, and DDX3X. In addition, the contribution of copy number variations (CNVs) to CVJ malformation was investigated and three pathogenic CNVs were identified in three affected individuals. To further dissect the complex mutational architecture of CVJ malformation, we performed a gene-based rare variant association analysis utilizing 4371 in-house exomes as control. Rare variants in LGI4 (carrier rate = 3.26%, p = 3.3 × 10-5) and BEST1 (carrier rate = 5.43%, p = 5.77 × 10-6) were identified to be associated with CVJ malformation. Furthermore, gene set analyses revealed that extracellular matrix- and RHO GTPase-associated biological pathways were found to be involved in the etiology of CVJ malformation. Overall, we comprehensively dissected the genetic underpinnings of CVJ malformation and identified several novel disease-associated genes and biological pathways.

A Hybrid Rule-Based and Data-Driven Approach to Illegal Transshipment Identification with Interpretable Behavior Features.

Illegal transshipment of maritime ships is usually closely related to illegal activities such as smuggling, human trafficking, piracy plunder, and illegal fishing. Intelligent identification of illegal transshipment has become an important technical means to ensure the safety of maritime transport. However, due to different geographical environments, legal policies and regulatory requirements in each sea area, there are differences in the movement characteristics and geographical distribution of illegal transshipment behavior in different time and space. Moreover, in areas with dense traffic flow, normal navigation behavior can easily be identified as illegal transshipment, resulting in a high rate of misidentification. This paper proposes a hybrid rule-based and data-driven approach to solve the problem of missing identification in fixed threshold methods and introduces a traffic density feature to reduce the misidentification rate in dense traffic areas. The method is both interpretable and adaptable through unsupervised clustering to get suitable threshold distribution combination for regulatory sea areas. The evaluation results in two different sea areas show that the proposed method is applicable. Compared with other widely used identification methods, this method identifies more illegal transshipment events, which are highly suspicious, and gives warning much earlier. The proposed method can even filter out misidentification events from compared methods' results, which account for more than half of the total number.

Functional characteristics of a broad spectrum of TBX6 variants in Mayer-Rokitansky-Küster-Hauser syndrome.

PURPOSE: Mayer-Rokitansky-Küster-Hauser syndrome (MRKHS) is characterized by congenital absence of the uterus, cervix, and the upper part of the vagina in females. Whole-gene deletion and loss-of-function variants in TBX6 have been identified in association with MRKHS. We aimed to expand the spectrum of TBX6 variants in MRKHS and explore the biological effect of the variant alleles. METHODS: Rare variants in TBX6 were called from a combined multiethnic cohort of 622 probands with MRKHS who underwent exome sequencing or genome sequencing. Multiple in vitro functional experiments were performed, including messenger RNA analysis, western blotting, transcriptional activity assay, and immunofluorescence staining. RESULTS: We identified 16 rare variants in TBX6 from the combined cohort, including 1 protein-truncating variant reported in our previous study and 15 variants with unknown effects. By comparing the prevalence of TBX6 variants in the Chinese MRKHS cohort vs 1038 female controls, we observed a significant mutational burden of TBX6 in affected individuals (P = .0004, odds ratio = 5.25), suggesting a causal role of TBX6 variants in MRKHS. Of the 15 variants with uncertain effects, 7 were shown to induce a loss-of-function effect through various mechanisms. The c.423G>A (p.Leu141=) and c.839+5G>A variants impaired the normal splicing of TBX6 messenger RNA, c.422T>C (p.Leu141Pro) and c.745G>A (p.Val249Met) led to decreased protein expression, c.10C>T (p.Pro4Ser) and c.400G>A (p.Glu134Lys) resulted in perturbed transcriptional activity, and c.356G>A (p.Arg119His) caused protein mislocalization. We observed incomplete penetrance and variable expressivity in families carrying deleterious variants, which indicates a more complex genetic mechanism than classical Mendelian inheritance. CONCLUSION: Our study expands the mutational spectrum of TBX6 in MRKHS and delineates the molecular pathogenesis of TBX6 variants, supporting the association between deleterious variants in TBX6 and MRKHS.

Lane-based Distance-Velocity model for evaluating pedestrian-vehicle interaction at non-signalized locations.

Pedestrian vehicle conflicts at non-signalized crosswalks are a world-wide safety concern. Although the "pedestrian priority" policy is applied in some regions to improve pedestrian safety, its effect needs further investigation. This study proposes the Lane-based Distance-Velocity model (LDV) to investigate pedestrian-vehicle interaction at non-signalized crosswalks. Compared with the DV model, the LDV model considers the lateral distance between vehicles and pedestrians. Therefore, the LDV model extends the application of the DV model by allowing it to be applied not only on one-lane streets to multi-lane streets. The conflict severities of pedestrian-vehicle interaction in the LDV model are classified into four categories: safe-passage, mild-interaction, potential-conflict and potential-collision. Based on that, pedestrian crossing decisions are graded as safe-crossing, risky-crossing, and dangerous-crossing. The experimental data are collected at a non-signalized crosswalk through drone footage collected in Xi'an City (China) with a Machine Vision Intelligent Algorithm. The model is tested through a case study to evaluate pedestrian crossing safety when interacting with private cars and taxis. Results from the case study suggest that the proposed model works well in the pedestrian-vehicle interaction analysis. Firstly, 87.9% of drivers are willing to provide right-of-way to pedestrians when they have enough time to react and yield. Then, both the DV model and LDV model have reached consistent conclusions: the deliberate violation rate (DVR) of taxi drivers is 22.64%, which is double that of private car drivers. Last, taxis commit a higher percentage of pedestrians' dangerous or risky crossing situations than private cars. Relevant government departments can utilize the results of this study to manage urban traffic better and improve pedestrian safety.

Identification of Novel FBN2 Variants in a Cohort of Congenital Contractural Arachnodactyly.

Congenital contractural arachnodactyly (CCA) is a rare autosomal dominant disorder of connective tissue characterized by crumpled ears, arachnodactyly, camptodactyly, large joint contracture, and kyphoscoliosis. The nature course of CCA has not been well-described. We aim to decipher the genetic and phenotypic spectrum of CCA. The cohort was enrolled in Beijing Jishuitan Hospital and Peking Union Medical College Hospital, Beijing, China, based on Deciphering disorders Involving Scoliosis and COmorbidities (DISCO) study (http://www.discostudy.org/). Exome sequencing was performed on patients' blood DNA. A recent published CCA scoring system was validated in our cohort. Seven novel variants and three previously reported FBN2 variants were identified through exome sequencing. Two variants outside of the neonatal region of FBN2 gene were found. The phenotypes were comparable between patients in our cohort and previous literature, with arachnodactyly, camptodactyly and large joints contractures found in almost all patients. All patients eligible for analysis were successfully classified into likely CCA based on the CCA scoring system. Furthermore, we found a double disease-causing heterozygous variant of FBN2 and ANKRD11 in a patient with blended phenotypes consisting of CCA and KBG syndrome. The identification of seven novel variants broadens the mutational and phenotypic spectrum of CCA and may provide implications for genetic counseling and clinical management.

Exome-wide Analysis of De Novo and Rare Genetic Variants in Patients With Brain Arteriovenous Malformation.

BACKGROUND AND OBJECTIVES: Brain arteriovenous malformation (bAVM) is a congenital disorder and a leading cause of hemorrhagic stroke. Germline genetic variants play an essential role in the pathogenesis of bAVM. However, the biological relevance of disease-associated genes identified in previous studies is elusive. In this study, we aim to systematically investigate the contribution of germline variants to bAVM and explore the critical molecular pathways underlying the pathogenesis of bAVM. METHODS: Probands with sporadic bAVM were consecutively recruited into this study from November 2015 to November 2018 and underwent exome sequencing. The controls were aggregated from individuals who were not known to have vascular malformation and underwent exome sequencing for clinical or research purposes. The retained control dataset included 4,609 individuals, including 251 individuals with parental samples sequenced. We first analyzed de novo variants in cases and controls and performed a pathway enrichment analysis. A gene-based rare variant association analysis was then performed to identify genes whose variants were significantly enriched in cases. RESULTS: We collected an exome-sequenced bAVM cohort consisting of 152 trios and 40 singletons. By first focusing on de novo variants, we observed a significant mutational burden of likely gene-disrupting variants in cases vs controls. By performing a pathway enrichment analysis of all nonsynonymous de novo variants identified in cases, we found the angiopoietin-like protein 8 (ANGPTL8) regulatory pathway to be significantly enriched in patients with bAVM. Through an exome-wide rare variant association analysis utilizing 4,394 in-house exome data as controls, we identified SLC19A3 as a disease-associated gene for bAVM. In addition, we found that the SLC19A3 variants in cases are preferably located at the N' side of the SLC19A3 protein. These findings implicate a phenotypic expansion of SLC19A3-related disorders with a domain-specific effect. DISCUSSION: This study provides insights into the biological basis of bAVM by identifying novel molecular pathways and candidate genes.

Expanding the mutation and phenotype spectrum of MYH3-associated skeletal disorders.

Pathogenic variants in MYH3 cause distal arthrogryposis type 2A and type 2B3 as well as contractures, pterygia and spondylocarpotarsal fusion syndromes types 1A and 1B. These disorders are ultra-rare and their natural course and phenotypic variability are not well described. In this study, we summarize the clinical features and genetic findings of 17 patients from 10 unrelated families with vertebral malformations caused by dominant or recessive pathogenic variants in MYH3. Twelve novel pathogenic variants in MYH3 (NM_002470.4) were identified: three of them were de novo or inherited in autosomal dominant way and nine were inherited in autosomal recessive way. The patients had vertebral segmentation anomalies accompanied with variable joint contractures, short stature and dysmorphic facial features. There was a significant phenotypic overlap between dominant and recessive MYH3-associated conditions regarding the degree of short stature as well as the number of vertebral fusions. All monoallelic variants caused significantly decreased SMAD3 phosphorylation, which is consistent with the previously proposed pathogenic mechanism of impaired canonical TGF-ß signaling. Most of the biallelic variants were predicted to be protein-truncating, while one missense variant c.4244T>G,p.(Leu1415Arg), which was inherited in an autosomal recessive way, was found to alter the phosphorylation level of p38, suggesting an inhibition of the non-canonical pathway of TGF-ß signaling. In conclusion, the identification of 12 novel pathogenic variants and overlapping phenotypes in 17 affected individuals from 10 unrelated families expands the mutation and phenotype spectrum of MYH3-associated skeletal disorders. We show that disturbances of canonical or non-canonical TGF-ß signaling pathways are involved in pathogenesis of MYH3-associated skeletal fusion (MASF) syndrome.

A genotype-first analysis in a cohort of Mullerian anomaly.

Müllerian anomaly (M.A.) is a group of congenital anatomic abnormalities caused by aberrations of the development process of the Müllerian duct. M.A. can either be isolated or be involved in Mendelian syndromes, such as Dandy-Walker syndrome, Holt-Oram syndrome and Bardet-Biedl syndrome, which are often associated with both uterus and kidney malformations. In this study, we applied a genotype-first approach to analyze the whole-exome sequencing data of 492 patients with M.A. Six potential pathogenic variants were found in five genes previously related to female urogenital deformities (PKD1, SON, SALL1, BMPR1B, ITGA8), which are partially overlapping with our patients' phenotypes. We further identified eight incidental findings in seven genes related to Mendelian syndromes without known association with reproductive anomalies (TEK, COL11A1, ANKRD11, LEMD3, DLG5, SPTB, BMP2), which represent potential phenotype expansions of these genes.

Aberrant interaction between mutated ADAMTSL2 and LTBP4 is associated with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is a complex spinal deformity with a prevalence of 1%-3%. Genetic factors have been associated with the etiology of AIS. However, previous studies mainly focused on common single nucleotide polymorphisms which confer modest disease risk. Recently, rare variants in FBN1 and other extracellular matrix genes have been implicated in AIS, suggesting a potential overlapping disease etiology between AIS and hereditary connective tissue disorders (HCTD). In this study, we systematically analyzed rare variants in a set of HCTD-related genes in 302 AIS patients who underwent exome sequencing. We firstly focused on pathogenic variants based on a monogenic inheritance and identified nine disease-associated variants in FBN1, COL11A1, COL11A2 and TGFBR2. We then explored the potential interactions between variants in different genes based on the case-control statistics. We identified three ADAMTSL2-LTBP4 variant pairs in three AIS patients and none in controls. Furthermore, we revealed that the variant pairs identified in these genes could affect the interaction between ADAMTSL2 and LTBP4 and upregulate TGF-ß signaling pathway in human fibroblasts. Our findings implicate that the aberrant interaction between mutated ADAMTSL2 and LTBP4 was associated with AIS.

The Mutational Landscape of PTK7 in Congenital Scoliosis and Adolescent Idiopathic Scoliosis.

Depletion of ptk7 is associated with both congenital scoliosis (CS) and adolescent idiopathic scoliosis (AIS) in zebrafish models. However, only one human variant of PTK7 has been reported previously in a patient with AIS. In this study, we systemically investigated the variant landscape of PTK7 in 583 patients with CS and 302 patients with AIS from the Deciphering Disorders Involving Scoliosis and COmorbidities (DISCO) study. We identified a total of four rare variants in CS and four variants in AIS, including one protein truncating variant (c.464_465delAC) in a patient with CS. We then explored the effects of these variants on protein expression and sub-cellular location. We confirmed that the c.464_465delAC variant causes loss-of-function (LoF) of PTK7. In addition, the c.353C>T and c.2290G>A variants identified in two patients with AIS led to reduced protein expression of PTK7 as compared to that of the wild type. In conclusion, LoF and hypomorphic variants are associated with CS and AIS, respectively.

Disruptive NADSYN1 Variants Implicated in Congenital Vertebral Malformations.

Genetic perturbations in nicotinamide adenine dinucleotide de novo (NAD) synthesis pathway predispose individuals to congenital birth defects. The NADSYN1 encodes the final enzyme in the de novo NAD synthesis pathway and, therefore, plays an important role in NAD metabolism and organ embryogenesis. Biallelic mutations in the NADSYN1 gene have been reported to be causative of congenital organ defects known as VCRL syndrome (Vertebral-Cardiac-Renal-Limb syndrome). Here, we analyzed the genetic variants in NADSYN1 in an exome-sequenced cohort consisting of patients with congenital vertebral malformations (CVMs). A total number of eight variants in NADSYN1, including two truncating variants and six missense variants, were identified in nine unrelated patients. All enrolled patients presented multiple organ defects, with the involvement of either the heart, kidney, limbs, or liver, as well as intraspinal deformities. An in vitro assay using COS-7 cells demonstrated either significantly reduced protein levels or disrupted enzymatic activity of the identified variants. Our findings demonstrated that functional variants in NADSYN1 were involved in the complex genetic etiology of CVMs and provided further evidence for the causative NADSYN1 variants in congenital NAD Deficiency Disorder.

Whole Exome Sequencing Uncovered the Genetic Architecture of Growth Hormone Deficiency Patients.

Purpose: Congenital growth hormone deficiency (GHD) is a rare and etiologically heterogeneous disease. We aim to screen disease-causing mutations of GHD in a relatively sizable cohort and discover underlying mechanisms via a candidate gene-based mutational burden analysis. Methods: We retrospectively analyzed 109 short stature patients associated with hormone deficiency. All patients were classified into two groups: Group I (n=45) with definitive GHD and Group II (n=64) with possible GHD. We analyzed correlation consistency between clinical criteria and molecular findings by whole exome sequencing (WES) in two groups. The patients without a molecular diagnosis (n=90) were compared with 942 in-house controls for the mutational burden of rare mutations in 259 genes biologically related with the GH axis. Results: In 19 patients with molecular diagnosis, we found 5 possible GHD patients received known molecular diagnosis associated with GHD (NF1 [c.2329T>A, c.7131C>G], GHRHR [c.731G>A], STAT5B [c.1102delC], HRAS [c.187_207dup]). By mutational burden analysis of predicted deleterious variants in 90 patients without molecular diagnosis, we found that POLR3A (p = 0.005), SUFU (p = 0.006), LHX3 (p = 0.021) and CREB3L4 (p = 0.040) represented top genes enriched in GHD patients. Conclusion: Our study revealed the discrepancies between the laboratory testing and molecular diagnosis of GHD. These differences should be considered when for an accurate diagnosis of GHD. We also identified four candidate genes that might be associated with GHD.

Novel FGFR1 Variants Are Associated with Congenital Scoliosis.

FGFR1 encodes a transmembrane cytokine receptor, which is involved in the early development of the human embryo and plays an important role in gastrulation, organ specification and patterning of various tissues. Pathogenic FGFR1 variants have been associated with Kallmann syndrome and hypogonadotropic hypogonadism. In our congenital scoliosis (CS) patient series of 424 sporadic CS patients under the framework of the Deciphering disorders Involving Scoliosis and COmorbidities (DISCO) study, we identified four unrelated patients harboring FGFR1 variants, including one frameshift and three missense variants. These variants were predicted to be deleterious by in silico prediction and conservation analysis. Signaling activities and expression levels of the mutated protein were evaluated in vitro and compared to that of the wild type (WT) FGFR1. As a result, the overall protein expressions of c.2334dupC, c.2339T>C and c.1261A>G were reduced to 43.9%, 63.4% and 77.4%, respectively. By the reporter gene assay, we observed significantly reduced activity for c.2334dupC, c.2339T>C and c.1261A>G, indicating the diminished FGFR1 signaling pathway. In conclusion, FGFR1 variants identified in our patients led to only mild disruption to protein function, caused milder skeletal and cardiac phenotypes than those reported previously.

Heterozygous Recurrent Mutations Inducing Dysfunction of ROR2 Gene in Patients With Short Stature.

PURPOSE: ROR2, a member of the ROR family, is essential for skeletal development as a receptor of Wnt5a. The present study aims to investigate the mutational spectrum of ROR2 in children with short stature and to identify the underlying molecular mechanisms. METHODS: We retrospectively analyzed clinical phenotype and whole-exome sequencing (WES) data of 426 patients with short stature through mutation screening of ROR2. We subsequently examined the changes in protein expression and subcellular location in ROR2 caused by the mutations. The mRNA expression of downstream signaling molecules of the Wnt5a-ROR2 pathway was also examined. RESULTS: We identified 12 mutations in ROR2 in 21 patients, including 10 missense, one nonsense, and one frameshift. Among all missense variants, four recurrent missense variants [c.1675G > A(p.Gly559Ser), c.2212C > T(p.Arg738Cys), c.1930G > A(p.Asp644Asn), c.2117G > A(p.Arg706Gln)] were analyzed by experiments in vitro. The c.1675G > A mutation significantly altered the expression and the cellular localization of the ROR2 protein. The c.1675G > A mutation also caused a significantly decreased expression of c-Jun. In contrast, other missense variants did not confer any disruptive effect on the biological functions of ROR2. CONCLUSION: We expanded the mutational spectrum of ROR2 in patients with short stature. Functional experiments potentially revealed a novel molecular mechanism that the c.1675G > A mutation in ROR2 might affect the expression of downstream Wnt5a-ROR2 pathway gene by disturbing the subcellular localization and expression of the protein.

Exome sequencing reveals genetic architecture in patients with isolated or syndromic short stature.

Short stature is among the most common endocrinological disease phenotypes of childhood and may occur as an isolated finding or in conjunction with other clinical manifestations. Although the diagnostic utility of clinical genetic testing in short stature has been implicated, the genetic architecture and the utility of genomic studies such as exome sequencing (ES) in a sizable cohort of patients with short stature have not been investigated systematically. In this study, we recruited 561 individuals with short stature from two centers in China during a 4-year period. We performed ES for all patients and available parents. All patients were retrospectively divided into two groups: an isolated short stature group (group I, n = 257) and an apparently syndromic short stature group (group II, n = 304). Causal variants were identified in 135 of 561 (24.1%) patients. In group I, 29 of 257 (11.3%) of the patients were solved by variants in 24 genes. In group II, 106 of 304 (34.9%) patients were solved by variants in 57 genes. Genes involved in fundamental cellular process played an important role in the genetic architecture of syndromic short stature. Distinct genetic architectures and pathophysiological processes underlie isolated and syndromic short stature.

Variants Affecting the C-Terminal of CSF1R Cause Congenital Vertebral Malformation Through a Gain-of-Function Mechanism.

CSF1R encodes the colony-stimulating factor 1 receptor which regulates the proliferation, differentiation, and biological activity of monocyte/macrophage lineages. Pathogenic variants in CSF1R could lead to autosomal dominant adult-onset leukoencephalopathy with axonal spheroids and pigmented glia or autosomal recessive skeletal dysplasia. In this study, we identified three heterozygous deleterious rare variants in CSF1R from a congenital vertebral malformation (CVM) cohort. All of the three variants are located within the carboxy-terminal region of CSF1R protein and could lead to an increased stability of the protein. Therefore, we established a zebrafish model overexpressing CSF1R. The zebrafish model exhibits CVM phenotypes such as hemivertebral and vertebral fusion. Furthermore, overexpression of the mutated CSF1R mRNA depleted of the carboxy-terminus led to a higher proportion of zebrafish with vertebral malformations than wild-type CSF1R mRNA did (p = 0.03452), implicating a gain-of-function effect of the C-terminal variant. In conclusion, variants affecting the C-terminal of CSF1R could cause CVM though a potential gain-of-function mechanism.

Whole-genome methylation analysis reveals novel epigenetic perturbations of congenital scoliosis.

Congenital scoliosis (CS) is a congenital disease caused by malformations of vertebrae. Recent studies demonstrated that DNA modification could contribute to the pathogenesis of disease. This study aims to identify epigenetic perturbations that may contribute to the pathogenesis of CS. Four CS patients with hemivertebra were enrolled and underwent spine correction operations. DNA was extracted from the hemivertebrae and spinal process collected from the specimen during the hemivertebra resection. Genome-wide DNA methylation profiling was examined at base-pair resolution using whole-genome bisulfite sequencing (WGBS). We identified 343 genes with hyper-differentially methylated regions (DMRs) and 222 genes with hypo-DMRs, respectively. These genes were enriched in the mitogen-activated protein kinase (MAPK) signaling pathway, calcium signaling pathway, and axon guidance in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and were enriched in positive regulation of cell morphogenesis involved in differentiation, regulation of cell morphogenesis involved in differentiation, and regulation of neuron projection development in Biological Process of Gene Ontology (GO-BP) terms. Hyper-DMR-related genes, including IGHG1, IGHM, IGHG3, RNF213, and GSE1, and hypo DMR-related genes, including SORCS2, COL5A1, GRID1, RGS3, and ROBO2, may contribute to the pathogenesis of hemivertebra. The aberrant DNA methylation may be associated with the formation of hemivertebra and congenital scoliosis.